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1.  Introduction 

It  is  known  that  a  two  dimensional  flow  of  inviscid  fluid  is 
unstable  with  respect  to  small  disturbances  if  the  velocity  profile  has 
an  inflexion  point  in  it.  Rayleigh  in  1880,  [1]  showed  that  the  con¬ 
dition  was  necessary,  and  Tollmien  [2]  showed  that  under  certain  con¬ 
ditions,  that  it  was  sufficient. 

The  purpose  of  this  report  is  to  investigate  the  effects  of  com¬ 
pressibility  on  the  stability  of  a  free  shear  layer,  whose  profiles 
have  points  of  inflexion  on  them.  The  hope  is  that  the  mechanisms 
involved  in  this  type  of  instability  will  explain  the  dissimilarites 
observed  between  the  stabilities  of  subsonic  and  supersonic  flows. 

As  a  starting  point,  we  use  Michalkes'  work  [3]  that  investigates 
the  stability  of  the  hyperbolic  tangent  velocity  profile  in  an  incom¬ 
pressible  flow.  We  use  this  profile  first  to  confirm  (for  our  own 
comfort)  Michalkes1  results  in  the  incompressible  care,  and  then  as  a 
vehicle  on  which  to  complete  the  compressible  study.  The  point  of  using 
this  profile  is  that  it  can  be  analytically  manipulated,  allowing  a 
relatively  simple  computational  code  to  be  used. 

Before  doing  the  compressible  study  proper,  we  first  eliminated  the 
possibilty  of  viscosity  being  a  significant  factor.  This  was  done  by 
demonstrating,  that,  for  the  incompressible  case,  viscosity  has  only  a 
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slight  attenuating  effect  on  the  type  of  instability  being  studied.  This 

3 

was  the  case  for  the  range  of  Reynolds  numbers  10  <  Re  < 


where  v(y)  is  the  imposed  velocity  profile,  is  the  free  stream  Mach 
number,  (J)(y)  is  the  complex  amplitude  of  the  disturbance  and  a  is  the 
complex  wave  number  of  the  disturbance.  A  derivation  is  sketched  in 
Appendix  1.  Boundary  conditions  on  (2.1)  are  derived  by  supposing  that 
for  an  unbounded  velocity  profile,  any  finite  disturbance  will  vanish  at 


infinity.  Thus,  for  boundary  conditions  on  (2.1),  we  have 

<J>(-°°)  =  -  o  (2.2). 

Equation  (2.1)  is  derived  by  assuming  a  disturbance  of  the  form 
A  -*(y)  -  ct) 

Thus,  the  disturbance  is  fully  described  when  4>  (y)  ,  ot  and  c  are  known. 


Clearly,  equation  (2.1)  with  (2.2)  is  insufficient  to  determine  all  of 
these  quantities,  and  so  some  further  restriction  is  necessary.  The 
classical  viewpoint,  [4]  is  that  such  disturbances  were  temporal  by 
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nature,  which  is  equivalent  to  saying  that  a  is  real  (i.e,  spatially 
bounded).  Thus,  for  a  specified  wave  number  a  (2.1)  with  (2.2)  reduces 
essentially  to  a  characteristic  value  problem,  with  c  as  the  characteristic 
value  and  <f(y)  as  the  characteristic  function.  In  recent  years,  however, 
the  view  has  formed  that  the  disturbances  are  spatial  in  character,  in 
which  case  ac  is  taken  to  be  a  real  quantity.  In  this  case  if  a  =  +  ja^ 

and  c  =  c_  +  jcT,  then  ac  -  real  implies  a  c  +  a  c  -  0  or 

K  1  K  1  IK 

°i  =  _aRcI/cR‘  (2'3) 

Thus,  for  a  specified  a  and  aT  given  by  (2.3),  and  (2.1)  with  (2.2) 

K  X 

again  becomes  a  more  or  less  conventional  characteristic  value/function 
problem  in  (c,  <J))  respectively. 

Both  approaches  were  considered  in  this  study,  and  the  conclusions 
concerning  the  behavior  of  the  instabilities,  with  varying  Mach  number, 
were  the  same  with  both  approaches.  In  the  light  of  this  coincidence  we 
chose  to  consider  the  disturbance  as  temporal,  mainly  to  facilitate  com¬ 
parison  with  published  work  (Michalke,  [3]). 

Thus,  the  mathematical  problem  as  so  far  defined  is  as  follows:  for 
a  specified  range  of  real  a's,  solve  the  characteristic  value  problem 
defined  at  (2,1)  and  (2.2).  Repeat  the  process  for  a  range  of  M^  (free 
stream  Mach  number) . 

3.  Organization  of  Equations  into  a  More  Convenient  Form. 

The  characteristic  value  problem  as  defined  at  (2.1)  and  (2.2)  is  not 
easily  solved,  but  a  solution  can  be  made  easier  by  following  the  procedure 
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Michalke  used  for  his  incompressible  study.  Putting  =  0  into  (2.1) 
gives,  after  rearrangement 


(v  -  c)(( J>"  -  0t2cJ) )  -  v"<J>  =  0  (3.1a) 

with 

<t>(“°°)  =  (+00)  =  0  (3.1b) 

This  is  Rayleighs’  stability  equation,  the  subject  of  Michalkes' 
study.  He  reduced  (3.1a)  from  a  second  order  equation  into  a  first  order 
equation  using  the  transformation 

<Ky)  *  exp  [  J 9dy ]  (3.2) 

to  obtain  from  (3.1a)  the  equation 

it 

0"  +  e2  -  a2  -  ( - - - )  =  0  (3.3) 

V  -  c 

To  obtain  the  boundary  condition(s)  on  0,  we  use  the  facts  that 

v"(_oo)  =  v"(+°°)  =  0  in  (3.1a).  Thus,  at  y  -  +  °°  (3.1a)  gives 

4>"  -  ot2cJ)  =  0  (3.4) 

yielding  at  y  =  +  00 . 

4> (y)  *  exp[+  ay]  (3.5) 

Since  we  know  that  as  y  *>  then  <J>(y)  0  then  (3.5)  gives  as  y  ■> 

then  4>(y)~*  exp[-  ay]  (assuming  a  >  0)  and  likewise,  as  y  +  f00  then 
(y)  =  exp(+  ay).  Comparison  of  these  two  solutions  with  (3.2)  yield 
that  as 

y  -►  -00,  then  0  ->  -Ha 

and  as  (3.6) 

y  -►  +°°,  then  0  ->  -a 

Applying  (6)  onto  the  equation  of  the  compressible  case  for  M  >  0 

and  following  similar  arguments  for  boundary  conditions  we  obtain  the  system 

(v  -  c)[g(0"  +  02)  -  g'0]  »  a2g2(v  -  c)(l  +  0.2M2)  +  v”  g  -  v'g"  (3.7) 
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with  the  corresponding  boundary  conditions 

0(_co)  s  aSign (a)  /gT-°°)  (3,8) 

Q(+oo)  =  -aSign(a)  /I  +  0.2M^  /g(+°°~) 

Here,  the  function  Sign  (a)  is  used  to  account  for  the  (non-physical) 
mathematical  possibility  of  a  being  negative.  See  Appendix  2  for  details. 

Superficially,  the  system  (3.7)  and  (3.8)  is  overdetermined,  in  that 
(3.7)  is  first  order,  and  yet  (3.8)  consists  of  two  boundary  conditions, 
to  be  satisfied  by  (3.2),  but  we  must  keep  in  mind  that  c  is  a  parameter 
at  our  disposal,  and  may  be  adjusted  to  ensure  the  second  (boundary) 
condition,  given  the  first. 

The  complex  system  of  (3.7)  and  (3.8)  yields  the  pair  of  simultaneous 


equations . 


=  - (02  -  02)  +  (AC  +  B°)  +  (AD  -  BC)0I  +  (AR1  -  BR2) 

R  1  - - — -  - 


(A2  +  B2) 


(A2  +  B2) 


<30-  (BC  -  AD)  0D  +  (BD  +  AC)0  r  (BR,  +  AR„) 

~  »  -20  6  +  - - - i  +  - i - — 

'  y  R  1  (A2  +  B2)  (A2  +  B2) 


where 


0  5  eR  +  jex 


(3.9) 


c  =  cR  +  jcx 
8  —  8r  j  B-j- 


A  =  gR(v  -  cR)  +  8i  Cj. 

B  =  -gj(v  -  cR)  +  8rCi 

C  =  8'r(v  "  CR>  +  8icI 


(3.10) 


"  'iii  .  r  ... 


d  =  -  gyv  -  cR)  +  g'RCl 

Rx  =  ct2(l  +  0. 2M2)  [  (gR  -  gj)(v  -  cR)  +  2clgRgl] 
+  v"gR  ~  v’g; 

R2  =  a2 (1  +  0.2M2) [-CjCg^  -  gj)  +  28r8x(v  “  CR)  J 

+  Vgj.  -  v'gj' 

with  the  boundary  conditions  as  y  -00 

0R  =  +a  Sign (a) (g2  +  g2)J  Cos  0/2 

6j  =  +a  Sign(a) (g2  +  g2)^  Sin  0/2 

A 

0  =  Principal  (Tan  1 gx/ gR) 

and 

as  y  +  +oo 

0R  =  -a  Sign(a)  (1  +  0.2M^j)'S  (g2  +  g2)^  Cos  0/2 


(3.10) 


0j  =  -a  Sign(a) (1  +  O^M2)4  (g2  +  gj)^  Sin  0/2  (3<1Q) 

/V 

0  =  Principal  (Tan  1 Sj/^r) * 

The  derivations  are  given  in  Appendix  2.  The  system  (3.9)  and  (3.10) 
is  still  not  in  a  suitable  form  to  be  solved  numerically,  because  the 
independent  variable  y  is  over  the  range  -°°  <  y  <  +°°.  To  eliminate  this 
problem,  we  introduce  the  independent  variable  transformation 

z  =  0.5(1  +  Tanh  (y)).  (3.11) 

This  transforms  the  range  <  y  <  +»  onto  0  <_  z  <_  1.  This  transformation 
has  a  simple  effect  on  (3.9),  viz:  if  we  label  the  two  right  hand  sides 
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(3.12) 


(3.9)  as  RHS^  and  RHS2  respectively,  then  (3.9)  becomes 
d9„  1 


=  -sr-Ti - TV  (RHS.) 

dz  2z(l  -  z)  1 


and 


d6  1 

dT  '  2z(1~~~zT  <RHV 


The  boundary  conditions  remain  the  same,  except  that  instead  oi 
being  evaluated  at  y  *  +  00 ,  they  are  evaluated  at  z  =  0,  1. 

It  is  clear  that  any  numerical  treatment  of  (3.12)  requires  care  at 
the  two  points  z  *  0  and  2*1.  This  treatment  is  described  in  Appendix  3. 

There  is  a  further  point  that  needs  to  be  clarified.  It  is  indicated 
in  Appendix  6  that  when  c =  0,  i.e.  at  the  neutral  points,  then  a  is 

necessarily  real  for  <_  2.5.  In  turn,  this  leads  to  a  0  so  that  in 

9  2  2  2 

(3.12),  the  denominator  A  +  B  reduces  to  g(v  -  c) .  Since  c  -  0.5  and 
0  <_  v  £  1,  then  clearly  the  region  about  v  =  c  needs  to  be  carefully 
treated.  This  treatment  is  given  in  Appendix  4. 


4.  Strategy 

Before  describing  in  detail  how  to  solve  the  system  (3.9)  and 

(3.10),  it  is  instructive  to  consider  the  objective.  We  are  primarily 

interested  in  the  change  in  behaviour  of  the  instabilities  associated 

with  a  given  velocity  profile  as  Mach  number  increases  in  the  range 

0  <_  <_  2.5,  (the  upper  limit  being  decided  by  considerations  to  be 

discussed  later).  From  [3]  it  can  be  seen  that  for  temporal  disturbances, 

for  a  given  wave  number  a,  the  growth  of  the  disturbance  is  determined 
aci 

by  e  .  The  intensity  of  a  disturbance  is  experimentally  observed  to 
depend  on  its  wave  length  A(»2Tr/a). 


ac 

It  is  instructive  to  plot  e  *  against  X  for  a  given,  fixed 
By  forming  a  set  of  such  curves  for  a  discrete  set  of  values  {M^}, 
a  visual  idea  may  be  obtained  of  the  effect  of  Mach  number  on  instability. 

Furthermore,  from  Appendix  1,  we  see  that  the  argument  of  the 
derivative  in  equation  (2.1)  is  directly  related,  in  a  fairly  simple 
fashion  to  the  complex  static  pressure  field  driving  the  disturbance. 

Since  static  pressure  measurements  may  be  (indirectly)  obtained  quite 
easily  from  experiment,  it  makes  sense  to  plot  the  pressure  fields  for 
various  Mach  numbers. 

5 .  Numerical  Procedure 

For  a  given  a  and  M,  the  system  (3.9)  and  (3.10)  has  the  basic  form 

9r  =  W  er  cr*  cr  z) 

®I  =  fI(0R’  0I»  CR’  ci’  Z) 

Vz  -  0)  =  eR0 

9I<Z  =  °>  =  eio 

where  the  pair  (cR,  c^)  has  to  be  such  that  the  following  conditions  are 
true: 

0R(z  -  1)  =  0R  (specified) 

1  (5.2) 

0T(z  =  1)  =  0  (specified) 

L1 

The  adopted  procedure  goes  as  follows:  For  a  given  (cD,  cT) ,  suppose  that 
the  solution  of  (5.1)  at  z  -  1  is  [0R(cR,  cj»  D*  6i^cr»  cj*  !)]•  Define 
the  functional  F  -  [0R(cR,  Cj,  1)  -  Q^]2  +  f0];(cR,  Cj,  1)  -  S^]2  (5.3) 
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The  original  characteristic  value  problem  now  reduces  to  the  problem 
of  minimizing  F  at  (5.3)  with  respect  to  cR  and  c^  Clearly,  -  0.0, 


and  this  state  occurs  when  (cR,  c^.)  is  such  that 

0O  (c  ,  c  ,  1)  *  0D 
K  K  I  Kj 


0I(CR’  CI*  *  9Ri 


(5.4) 


i.e.,  F  is  minimized  when  the  original  characteristic  value  problem  is 

solved.  Thus,  for  a  given  fixed  value  of  the  pair  (M^,  a)  we  can  generate 

the  corresponding  pair  (c  ,  cT).  To  generate  the  trajectory  (a,  c  ,  cT) 

for  a  fixed  M^,  we  advance  a  to  a  +  Aa  and  then  calculate  (cR  +  Ac^, 

Cj  +  A Cj).  The  process  is  repeated  until  a  sufficient  range  of  a  is 

covered.  As  explained  in  the  last  section,  we  are  interested  in  a  plot 

of  a  against  c^,  and  so  essentially,  we  form  a  projection  of  the  trajectory 

(a,  c  ,  c  )  onto  the  (a,  c  )  plane.  In  practice,  it  turns  out  that  c 
K  1  J.  K 

varies  only  very  slightly  along  such  a  trajectory  (0.6%  for  M  =  1), 

and  so  little  if  any  information  is  lost  by  this  particular  representation. 


5 . 1  Difficulties  Involved 

As  with  many  apparently  straight  forward  problems,  hidden  difficulties 
are  often  encountered.  The  present  case  is  no  exception.  The  whole 
method  hinged  upon  our  knowing  for  each  value  of  M  ,  just  one  point  on 
the  trajectory  (a,  cR,  c^) .  It  transpired  that  it  was  very  difficult  to 
locate  any  point  on  it.  The  solution  was  to  adopt  the  old  African  adage: 
"slowly,  slowly  catchee  monkee...",  and  to  proceed  as  follows.  If  the 
velocity  profile  is  assumed  to  be  the  hyperbolic  tangent  profile,  then 
at  «  0,  it  is  known  (Michalke)  that 

(a,  cR,  cx)  *  (1,  0.5,  0.0)  (5.5) 
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is  a  point  on  the  trajectory.  It  is  the  neutral  point  on  the  trajectory. 

We  next  assume  that  all  trajectories  corresponding  to  in  the  range 
0  <_  <_  2.5  have  such  neutral  points,  (a,  c^,  0).  Starting  from  =  0 

and  by  incrementing  by  some  AM^,  that  is  suitably  small,  a  trajectory 

•k 

(M  ,  a,  c  )  may  be  generated  for  c  =  0.  From  this  trajectory,  we  can 
00  K  I 

pick  the  neutral  point  (a,  c  ,  0)  for  any  of  the  trajectories  M  =  constant. 

i\  00 

Knowledge  of  this  point  then  allows  the  generation  of  the  trajectories 

(a,  c  ,  c  )  for  the  given  M  =  constant. 

K  1  00 

If  the  velocity  profile  is  not  the  hyperbolic  profile,  it  is  just 
assumed  that  the  neutral  point  for  =  0  lies  in  the  region  of 
(a,  cR,  Cj)  *  (1,  0.5,  0)  and  a  search  is  conducted.  Having  found  that  one 
point,  everything  proceeds  as  described  above. 

6.  Results 

The  hyperbolic  tangent  profile  is  given  by 

v(y)  =  0.5(1  +  Tanh  (y)).  (6.1) 

In  the  range  -°°  <  y  <  °°,  then  0  v(y)  <  +1.  At  y  =  0,  there  is  a  point 
of  inflexion.  The  curve  is  sketched  in  Figure  1. 

6.1  Effectiveness  of  Current  Procedure 

This  profile  was  chosen  because  a  study  on  the  instabilities  in  an 
incompressible  flow,  using  this  profile,  has  been  conducted  by  Michalke  [3]. 
Thus,  in  the  one  case  =  0,  his  results  form  a  ready  basis  for  comparison 
with  our  results,  giving  us  a  measure  of  the  effectiveness  of  the  presently 
used  procedure. 

*Note;  It  transpires  that  for  c^  *  0,  then  an  algebraic  relation 
A(Moo,  c^)  *  0  exists.  Details  are  explained  in  Appendix  4. 
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Table  1  displays  the  results  of  Michalke  and  those  of  the  present 
method  for  =  0.  It  is  readily  seen  that  (to  within  rounding  error) 
perfect  agreement  exists.  We  were  thus  able  to  proceed  with  a  certain 
degree  of  confidence  in  both  the  method  and  the  code. 

6 . 2  Effect  of  Viscosity 

The  next  step  was  to  determine  the  effect  of  Reynolds  number  on  the 
instabilities  in  the  flow  [5].  The  stability  equation  for  the  case  of  finite 
Reynolds  numbers,  and  =  0  is 

(v  -  c)(4)M  -  a240  +  -  2a24>”  +  a44>)  -  vM(J)  =  0  (b.2) 

aRe 

Application  of  the  transformation  procedure  described  in  Section  3  yields 
the  system 

A"  +  A2  -  a2  =  0 

^L(0"  -  a20)  +  (v  -  c)0  =  v"  (6.3) 

=  exp  (I  Ady) 

0(-oo)  =  0  (+°°)  =  0 

The  details  of  this  derivation  are  given  in  Appendix  5.  The  system  (23) 
was  solved  using  a  modified  version  of  that  procedure  described  in  Section 
5.  See  Appendix  5.  Results  are  displayed  only  for  the  worst  case, 

3 

namely  Re  *  10  ,  in  Table  2.  Column  (a)  gives  the  Re  =  00  case,  and  column 

3 

(b)  gives  the  Re  -  10  case.  It  is  plain  that,  even  at  this  relatively 
low  Reynolds  number,  the  viscous  attenuation  of  the  instabilities  is 
marginal . 

Consequently,  it  was  felt  that  viscosity  could  safely  be  ignored  as 
a  significant  mechanism  in  the  type  of  instabilities  currently  being 
considered. 
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6. 3  Effect  of  Compressibility 


Computations  were  performed  for  Mach  numbers  in  the  range  0  <_  <  *2.5. 

The  upper  was  forced  upon  us  because  at  =  2.3,  the  neutral  point 
trajectory  (a,  c  ,  0)  enters  the  complex  plane,  indicating  the  limit 
of  the  instability  under  consideration.  See  Appendix  6. 

6.3.1  Behavior  of  the  Neutral  Point 

Figure  2  displays  the  behavior  of  the  neutral  point  (a,  c^)  =  (a,  0) 
as  increases  in  the  range  0  <_  <  2.5.  Table  3  displays  the  same 

results  numerically. 

In  essence,  the  neutral  point  appears  to  move  monotonically  from  1.0 
to  0.0  as  moves  from  0.0  to  2.5,  indicating  a  progressive  decrease  in 
the  range  of  unstable  wavelengths  X(=  2u/a)  for  disturbances,  as  increases. 
This  last  statement  is  qualified,  because  there  seems  to  be  an  anomolous 
behavior  present  in  a  at  about  =  1.375,  not  shown  graphically,  indicat¬ 
ing  the  possible  presence  of  a  singularity  in  a.  Unfortunately,  the 
presently  used  code  did  not  behave  very  well  in  the  region  concerned,  and 
so  we  were  unable  to  readily  describe  the  events  occurring. 

From  the  gross  behavior  of  the  neutral  point,  two  conclusions  may  be 
drawn  -  the  results  indicate 

(a)  Only  large  wavelengths  are  unstable  (X  >  2tt6)  or  only  low 

frequencies  (f  <  c  / 2tt6 )  are  unstable. 

R 

(b)  At  high  Mach  numbers,  the  range  of  unstable  wavelengths  is  much 
less  that  at  low  Mach  numbers. 

These  results  appear  to  be  consistent  with  the  available  experimental 
data  [6,  7]. 


12 


6.3.2  Behavior  of  Instabilities  Over  All  Wavelengths 


These  results  are  displayed  over  the  Mach  number  values  =  0.0, 

0.5,  1.0  and  2.0.  Tables  4a,  4b,  4c,  and  4d  give  the  numerical  results 
for  the  four  respective  Mach  numbers.  Figures  3  gives  the  graphical 
displays.  There  are  three  conclusions  to  be  drawn: 

(a)  There  is  apparently  no  high  wavelength  cut-off  to  the  instability 
of  disturbances. 

(b)  The  magnitude  of  amplification  decreases  as  Mach  number  increases 

(c)  The  wavelength  at  which  maximum  amplication  occurs  increases  as 
Mach  number  increases. 

6.4  Determination  of  the  Eigen  Function 

Solution  of  the  characteristic  equation  produces  eigen  values  for  the 
complex  velocity,  c(a,  M^)  and  also  the  eigen  function,  <f> (a ,  Mm)  . 

The  function  of  most  interest  is  the  pressure  P(y)  since  comparison 
with  experimental  pressure  data  may  be  possible.  Referring  to  equation 
A. 1.9  in  Appendix  1 

yP  2  r  =  “  og  KV  “  ■  v  ^  <6-4) 

The  real  and  imaginary  values  of  V  are  shown  in  Figures  4  and  5  respectively. 
For  a  given  M^,  the  4>  plotted  is  the  one  whose  associated  c  gives  maximum 
amplification.  For  comparison  with  experiment  the  pressure  may  be  obtained 
as  follows: 

P'  *  Real  P(y)  exp[ja(x  -  ct)] 

The  two  functions  Real  (T)  and  Imaginary  (T) ,  (referred  to  as  R(T)  and 
1(f)  hereafter)  are  plotted  in  Figures  4  and  5  respectively. 


They  are  plotted  on  the  z-axis,  0  z  <_  1,  corresponding  to  a  y-axis  plot 

<  y  <  +°°.  Now,  since  as  y  -►  +«>  ,  <f>"(y)  ->  0  and  <J>(y)  -►  0,  then  from 

equation  (6.4),  R(T)  0  and  I(r)->  0  as  y  i.e.,  as  z  -►  0  or  z  1 , 

then  I (T)  0  and  R(T)  ->  0.  From  Figure  5,  1(0  is  seen  to  behave  correctly. 

From  Figure  4,  for  *  0.0  and  1.0,  R(0  also  is  seen  to  behave  correctly, 
but  for  M  =1.5  and  2.0  R(T)  appears  to  be  tending  to  some  other  value, 

<x> 

However,  at  z  =  0  and  z  =  1,  then  <p(z)  =  0  and  by  virture  of  (A./. 4), 

<p'(z)  -  0  at  z  =  0,  1.  It  follows  that  at  z  =  0,  1,  then  R(0  must  neces¬ 

sarily  also  be  zero.  This  implies  strange  behavior  for  R(T)  when  =  1.5 
and  M  =2.0,  which  quite  possibly  is  numerical  in  origin.  It  may  also  be 
connected,  somehow,  to  the  appearance  of  instabilities. 

7 .  Summary 

The  stability  of  a  compressible  free  shear  layer  has  been  investigated 
by  utilizing  the  linearized  equation  resulting  from  a  small  perturbation 
analysis  (Rayleigh  equation) .  This  eigen-value  problem  was  solved  numerically 
for  various  wave  numbers  (a)  and  Mach  number  for  a  mean  velocity  profile  of 
the  form  of  a  hyperbolic  tangent.  This  type  of  flow  (with  an  inflection 
point)  was  shown  to  be  unstable  for  low  values  of  wave  number  and  for  Mach 
numbers  below  2.5.  These  results  are  believed  to  be  useful  for  analyzing 
aerodynamics  instabilities  [8]  encountered  in  separated  flows  (buzz,  cavity 
resonance,  buffet,  etc,). 
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APPENDIX  A 


DERIVATION  OF  THE  GOVERNING  EQUATIONS 

An  inviscid,  compressible  unsteady  flow  is  assumed.  Let  u  represent 
the  streamwise  cartesian  velocity  component,  and  v  the  cross  flow  component, 
p,  p  and  T  are  the  static  pressure,  density  and  temperature  respectively. 

The  basic  equations  are:  Continuity,  two  momentum  and  energy. 


DP 

pDt 


+  V  •  V 


0 


Du  _  _  dp 


Dt 


3x 


D v  _  3p 


(A. 1.1) 


Dt 


3x 


Dp  ^  yDp 
pDt  pDt 


D 


where  —  »  substantial  derivative 

Eliminating  the  density  from  the  continuity  and  energy  equation 

produces  the  following: 

3u  3u  3u  1  3p  n 
9F+u^  +  v97+  p  3^  =  0 

3v  _I_  1  3p  r\  /a  i  o \ 

3F+u3^  +  vW  +  ?3y  (A-1’2) 

i£  +  ui£  +  v92.  +  Yo,1h.  +  =  o 

3t  u  3x  v  3y  Yp(3x  +  3y'  U 
Now  assume  that  the  flow  can  be  represented  as  a  disturbance  from  a 
steady  state  shear  flow  [4,  5]  in  which 
u(y)  *  streamwise  velocity 


v  «  0  *  cross  flow  velocity 
p  ■  static  pressure  -  constant 


(A. 1.3) 
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It  follows  that  -  77^"  -  0. 

dx  dx  dy 

Let  the  time  dependent  disturbed  flow  be  given  by 
u  =  u  +  u1 


v  =  v 


P  =  P  +  P? 


With  this  flow,  the  equations  (A. 1.2)  give 


9uf 

8t 


+  u 


9uf 

9x 


I 

p  9x 


0 


9vr 

9t 


+  u 


9v' 

9x 


+ 


1  3£* 

p 


0 


3r' 

at 


+  yp( 


-,9u* 


8x 


+  £')  =  0 
3y 


(A. 1.4) 


(A. 1.5) 


We  now  assume  periodic  disturbances  of  the  form 
u'  =  u(y)  exp[ja(x  -  ct) ] 

v'  =  v(y)  expt ja(x  -  ct)]  (A. 1.6) 

p’  =  p (y)  exp[jot(x  -  ct)] 

where  u,  v  and  p  are  assumed  to  be  complex.  (A. 1.6)  substituted  into 
(A. 1.5)  yields,  in  turn, 

A. 

ja(u  -  c)  u  +  u^v  *  -  ^  jot  (A.  1.7a) 

ja(u  -  c)v  =  ~  ~  Py  (A. 1.7b) 

ja(u  -  c)p  =  -  yp(jau  +  v  )  (A. 1.7c) 


We  seek  to  eliminate  p  and  u  from  these  three  equations.  From  (A. 1.7c) 

A 

jau  =  -  ja(u  -  c)  -§  -  vy  (A. 1.8) 

In  (A. 1.7a),  (A. 1.8)  gives 

(u  -  c)[-ja(u  -  c)  yjr  -  0y]  +  Uy$  “  -  i  ja 

1  /"  N2 

u  Q  -  v  (u  -  c)  *  -  jct[=- - 1  P 

y  7  w  YP 
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Use  p  =  p/RT,  and  we  get 


V>  “  °y(u  "  C)  .  p 

~T -  _ - 2  "  "  JOt  vp 

yRT  -  (u  -  c)Z  YP 


(A. 1.9) 


Now  eliminate  p  between  (A. 1.7b)  and  (A. 1.9)  and  we  get 


u  v  -  (u  -  c) v 

y _ 2. 

yRT  -  (u  -  c)4 


(ja) (YRf)  (u  ' 


a  y 


or 


u  v 

Jl  . 


-  (u  -  c)v 


yRT  -  (u  -  c)' 


a^u  - 


c)v 


yRT 


-J  y 


(A. 1.10) 


Now  use . 


yRT 


2  -2 

1  -  0.2  M  u 

o 

M2 


(A. 1.11) 


where 


M2  = 
o 


M 


1  +  0.2  M 


(A. 1.12) 


and  (A. 1.10)  becomes 


u  v  -  (u  -  c ) v 

_  y___  y 


2-2  2  -  2 

(1  -  0.2  M  u  )-M  (u-c) 

o  o 


-  q2(u  -  c) v 

(1-0.2  h2  n2) 


(A. 1.13) 


Putting 


M2  = 


m2  ~2 

M„  u 


1  +  0.2  M2  (1  -  u2) 


(A. 1.14) 


yields,  after  some  algebra 


u  v  -  (u  -  c)v 

JL  y 


(1  +  0.2  M2)-1  -  M2(G  -  c)2 


a2(l  +  0.2M2)(u  -  c)v  (A. 1.15) 


a  y 


17 


A  change  of  notation 

V  =  0 

U  =  v 

yields  equations  (2.1) 
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APPENDIX  B 


THE  DERIVATION  OF  BOUNDARY  CONDITIONS  FOR 
THE  TRANSFORMED  STABILITY  EQUATION 


The  un transformed  stability  equation  for  an  inviscid,  compressible  flow 


is  given  by 


fL  CH.'.  -  v,<t>.)  =  a2d  +  0.2M2)  (v  -  c)<}>  (A. 2.1) 


We  consider  this  equation  for  very  large  value  of  |y|  <  °°. 


As  y  then  v  0,  M  0,  g  1  -  M^c  and  (A.  2.1)  may  be  approximated 


i.e.  by 


-c4)n  =  a  (1  -  M  c  )(-c)4) 
o 


<p"  -  a2(l  -  M2c2)0  =  0 
o 


(A. 2. 2) 


(A. 2. 3) 


i.e.  as  y  -*  then  <P  -►  exp(+  a(l  -  M  c  )  2  y) 

—  o 


(A. 2. 4) 


Since  the  boundary  condition  on  $  at  is  $(-<»)  *  0,  then  (A.  2. 4)  must  be 


where 


(p  =  exp(aa(l  -  M^c^)^  y) 
a  =  Sign(Real(a(l  -  M2c2)^)) 


Since  by  (3.2) 


(A. 2. 5) 


(A. 2. 6) 


<p  =  exp(/0dy) 


(A. 2. 7) 


Then  comparison  of  (A. 2. 7)  with  (A. 2. 5)  yields  that  a  y  -+  then 


0  -*•  a  a(l  -  hTc  ) 2 


2  2  5s. 

A  „  -  \  2\ 


a  =  Sign(Real(a(l  -  JTc‘)a)) 


(A. 2. 8) 


80=  8R  +  j8j  =  (1  -  Mqc  ), 

8o  =  <gR  +  8I)iS  exPOTan^1(gj/8R) ) 


then. 


(A. 2. 9) 


Put  0  =  Tan  hg-r/go).  then  (A.2.8)  may  be  written 
1  i\ 

0  aa(g^  +  g^)4  exp ( j0/2) 
a  =  Sign(a  Cos  0/2) 


(A. 2. 10) 


Splitting  (A. 2. 10)  into  its  real  and  imaginary  parts  yields  the  boundary 
conditions  at 


eR  =  ctS(g2  +  g2 )h  Cos(0/2) 
61  =  aa(gR  +  g2)U  Sin (0/2) 


(A. 2. 11) 


a  =  Sign(a  Cos(0/2)) 

9  =  Tan_1(gI/gR) 

As  y  4  +-°°,  then  v  -*■  1,  M  -*■  and  g  -*■  g^  =  (1  +  0.2  M2)  ^  -  M2(l  -  c)2. 


<J>"  -  a2  g  (1  +  0.2  M2)<f>  =  0 


so  that  for  large  positive  y,  (A. 2.1)  may  be  approximated  by 

(A. 2. 12) 

Following  similar  arguments  to  those  used  to  deduce  (A. 2. 11),  we  arrive 
at  the  boundary  condition  at 


y  =  +°° 

0R  =  -  cca(g2  +  g2)^  Cos(0/2) 

0r  =  -  aa(g2  +  g2)^  Sin (0/2) 

A  /\ 

a  =  Sign(a  Cos(0/2)) 

6  =  Tan"1(gI/gr) 


(A. 2. 13) 


In  evaluating  0  *  Tan  (gj/g^),  care  must  be  taken  to  ensure  that  the 
correct  quadrant  Is  used. 
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APPENDIX  C 


Thus,  evaluated  at  z  =  0,  (A. 3. 4)  gives 

(i  +  9„)  eR  -  9^  -  1/2  fR 

9I®8  +  (1  +  eR)  -  1/2  f  j 


(A. 3.1) 


(A. 3. 2) 


(A. 3. 3) 


(A. 3. 4) 


)A.3.5) 
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Solving  (A. 3.5)  for  0  and  0  gives 

K  1 

(1  +  9r)  K  +  9I  *i 

e  (Z  =  o)  =  1/2 ( - -  —  1  -1-  > 

(i  +  e  r  +  ej 


z  =  0 


(i  +  e  )  f  -  e  f 

K  (Z  =  0)  =  1/2  (— - —  -  —■  — ) 

(i  +  e  r  + 


z  =  0 


Following  the  same  arguments  at  z  =  1  gives 


0R  (z  =  1)  =  —  1  / 2 (• 


(2  -  6R)  fI  -  6I  C, 

«  -  0R)2  +  e2 


Z  =  1 


(1  -  0  )  f  -  0  t 

0J  (2  -  l)  -  -1/2C - s  _£  -  1 — S) 


a  -  eE)2  ♦  e2 


Z  =  1 


(A . 3 . 6) 


(A. 3. 7) 
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Now  use 


(A. 4b. 10) 


2(2v  -  1)(1  +  0 . 2M  )  -  0 . 4  MM1  =  Q 


M2  = 


2  2 
M  v 


1  +  0.2M2(1  -  v2) 


(A. 4b. 11) 


and  we  get,  with  v  =  c 
M  +  0.2  M  c  M  =  (• 


o.i  m: 

2c  -  P  ~  Om  ~  0.8c^(c  -  1)M2) 


where 


M  *  1  +  0.2  (1  -  c2) 


Now  put  a  =  0.2  and  rearrange,  these  results 

-c  +  4c  -  2)  +  (2c  -  1)  =  0 


2  2 

o  ("C  +  2c  -  1)  +o(-  2 


This  yields 

a  =  _(2c  -  1) 


(A. 4b. 12) 


(A. 4b. 13) 


(A. 4b.  14) 


(c  -  1)‘ 


M  =  {1L2  C  -  H{> 

(c  -  ir 

Alternatively,  we  can  rearrange  this  to  give 

c  -CL*  0-2  m£)  -  (l  +  Q.2  m£)  [ 

0.2  M2 


(A. 4b. 15) 


(A. 4b. 16) 


The  implication  of  this  result  is  significant  namely,  that  for  a  given 

profile  v(y) ,  then  the  propagation  velocity  of  the  neutral  aistrubance 
uniquely  defined  by  the  Mach  number  at  infinity. 


IS 


APPENDIX  E 


MATHEMATICAL  AND  NUMERICAL  TREATMENT  OF  THE 
STABILITY  EQUATION  FOR  INCOMPRESSIBLE,  VISCOUS  FLOWS 


The  equation  is  given  by 

(v  -  c)(<j>"  -  a2<f))  +  j/aRe(<}>  -  2a%  "  +  a^<f>)  -  v"(J)  =  0 

Using  the  D-operator  rotation  this  may  be  written  as 

(v  -  c)(D2  -  a2)<p  +  j/aRe(D2  -  a2)  (D2  -  a2)<p  -  v"4>  =  0 

Put  4>  =  exp(jAdy)  in  A.  5.2,  and  we  get 

(v  -  c)(A’  +  A2  -  a2)  +  j/aR  (D2  -  a2) (A’  +  A2  -  a2)  -  v" 


Put 


A’  +  A2  -  a2  =  0 


and  the  equations  are  now 

2  2  2 
r  +  r  -  a  =  0 


j/aR  (0M  -  a^O)  +  (v  -  c)0 


Using  the  arguments  of  Section  3,  the  boundary  conditions  of  A  are 
A(-°°)  =  +ot,  A(+°°)  =  -a 
Using  these  in  (A. 5.4)  gives 
0(-°°)  *  0(-M»)  -  0 

The  total  system  it 

A’  +  A2  -  a2  -  0 

j/aRe(0"-  a20)  +  (v  -  c)0  -  v" 

0(_OO)  -  0 (+<*>)  -  o 

A(-oo)  »  -Kx 

where  c  is  such  that  A(+°°)  »  -a 
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(A. 5.1) 

(A. 5.2) 

=  0  (A. 5. 3) 

(A. 5. 4) 

(A. 5. 5) 

(A. 5. 6) 

(A. 5. 7) 

(A.5.8) 


The  numerical  procedure  used  was  essentially  that  described  in 
Section  5,  except  that  for  each  value  of  c,  the  differential  equation 
for  0  must  be  solved.  Note  that  as  the  differential  equation 

in  0  becomes  an  algebraic  relationship  in  0*  and  the  system  reverts  to 
the  inviscid  form. 
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APPENDIX  F 


THE  UPPER  LIMIT  ON  MACH  NUMBER 


In  this  appendix,  it  is  shown  that  a  necessary  condition  for  the 

neutral  point  (a,  c  )  =  (a,  0)  to  have  a  real  ax  is  that  M  <  2,5.  It 
I  00  — 

is  not  shown  that  it  is  also  sufficient,  because  the  numerical  calculations 


demonstrate  that  it  is. 

The  stability  equation  is  given  by 

d_((y  -  c)<fr'  -  v'4>)  =  a2(1  +  2m2)(v  _  c)(J) 


g  =  (1  +  0.2M2)  -  M2(v  -  c)2 

M2  =  M2v2/(1  +  0.2M2  (1  -  v2)) 


M2  =  M2/(l  +  0.2M2) 

O  co  00 7 


I 


(A. 6.1) 


Suppose  that  for  some  y  then  g  is  such  that  g(y)  <_  0  for  all 

*  2 
-c°  £  y  <  y  •  We  will  show  that  in  this  case  that  a  is  necessarily  complex. 

It  is  known  that  as  y  -►  then  g  constant,  -w  say.  Thus,  for 
*  2 

y  <<  y  then  g  -  -  w  <0.  Likewise,  v  -  0  and  vf  =0.  Hence,  (A. 6.1)  may 
be  written  as 

xtt  2 

<p  m  -w  a  4> 


for 

* 

“°°  <  y  <<:  y 

(A. 6. 2) 

and 

<K-°°)  =  o 

If  we  insist  on  non-trivial  solutions  to  (A. 6. 2),  then  a2  must  be 
complex,  otherwise,  the  equation  in  (A. 6. 2)  will  only  admit  sinusoidal 
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solutions  that  cannot  satisfy  the  boundary  condition  $(-°°)  =  0. 

2 

Thus,  to  obtain  a  real  a  ,  then  g  must  remain  positive  at  least  at 
its  end  points.  It  is  shown  in  Fig  j  that  g  is  wholly  positive  for 
0  <  Mqo  <  2*^»  that  S (~°°)  “  g(+cc)  =  0  at  Moo  =  2,5  and  that  g(”°°)  < 
g(+°°)  <  0  for  Moo  >  2.5. 

Thus,  M^  =  2.5  represents  an  upper  bound  on  Mach  numbers  that  allow 
real  values  of  a  at  the  neutral  points. 


APPENDIX  G 


USEFUL  RELATIONSHIPS  BETWEEN  $  AND  9,  FOR  COMPUTATIONAL  PURPOSES 

Equation  (3.2)  relates  0(y)  and  <(> (y )  using  <Ky)  =  exp[j6dy]  (A. 7.1) 

where 

-OO  <  y  <  ao  . 

For  computational  purposes  it  is  necessary  to  "fix"  a  bottom  limit  to  the 
integral  in  (A. 7.1).  Thus  more  precisely, 

y 

(y)  s  exp[  4>(s)ds]  (A. 7.2) 

o 

Fixing  the  bottom  limit  at  s  -  0  ensures  that  <f>(0)  =  1  i.e.,  <J>(y)  is  normal¬ 
ized  to  unity  at  the  midpoint  of  the  range  of  y. 

From  (A. 7. 2)  there  follows 

y  y 

<f>R(y)  =  exp[  0R  ds]  cos[ 

o  o 

y  y 

<t’I(y)  =  expt  eR  ds]  sin[ 


ex  ds] 


0J  ds] 


(A. 7. 3) 


dy 


9r  *r  -  9I  *1 


dy 


eR  *1  -  9I 


(A. 7. 4) 


In  the  z-plane,  (A. 7. 2)  becomes 

4>  ( z )  ■  exp [  - — — - —  ds]  (A. 7. 5) 

2s(l  -  s) 

0.5 

Relationships  corresponding  to  (A. 7. 3)  follow  directly,  and  hence,  values 

for  and  <p'T  can  be  found  on  the  z-plane  using  (A. 7.4). 

K  X 
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TABLE  1 


Comparison  of  Michalkes  Results  With  Those  of  the 
Present  Method,  for  Incompressible  Inviscid  Flow 


C1 

a(MICHALKE)  1 

1 

a (PRESENT  METHOD) 

0.5000 

0.0 

i 

0.0  1 

i 

0.3487 

0.2000 

0.2002 

i 

0.2133 

| 

0.4446  ; 

0.4443  j 

0.1442 

I 

0.6000 

1 

0.5995  1 

| 

0.0674 

0.8000 

j 

1  0.7995  1 

|  j 

0.0000 

1.0000 

0.9997  J 

TABLE  2 

Effect  of  Reynolds  Number  on  Amplification 
Of  Disturbances  for  Incompressible  Flow 


(a) 


a,  Re  -  00 


0.5000 

0.0 

0 

0.3487 

0.2002 

0.2016 

0.2133 

0.4443 

i  0.4404 

0.1442 

0.5995 

0.5939 

0.0674 

0.7995 

0.7920 

0.0000 

0.9997 

Not  computed 

TABLE  3 


Behavior  of  the  Neutral  Point 

(a,  C  ,  C  )  =  (a,  C  ,  0)  as  Mach  Number  Increases 
K  I  K 


M 

00 

a 

CR 

0.0 

0.9973 

0.500 

0.2 

0.9934 

0.501 

0.4 

0.9747 

0.504 

0.6 

0.9440 

0.509 

0.8 

0.9022 

0.515 

1.0 

0.8502 

0.523 

1.2 

0.7888 

0.532 

1.4 

0.7187 

0.541 

1.6 

0.6396 

0.552 

1.8 

0.5509 

0.562 

2.0 

0.4483 

0.573 

2.2 

0.3194 

0.584 

2.4 

0.0796 

0.505 

2.5 

Not  Computed 

0.600 
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TABLE  4a 

Amplication  Against  Wavelength  for  Mach  Number  0.0 


o 

II 

s8 

.0 

exp (aCj) 

X  ~  2tt/oi 

1.0000 

6.28 

1.0605 

8.21 

1.0701 

9.10 

1.0926 

10.88 

1.0947 

11.33 

1.0965 

11.82 

1.0978 

12.35 

1.0987 

12.90 

1.0992 

13.50 

1.0994 

14.14 

1.0992 

14.83 

1.0987 

15.57 

1.0978 

16.38 

1.0967 

17.25 

1.0952 

18.19 

1.0935 

!  19.22 

1.0914 

20.34 

1.0892 

21.58 

1.0866 

22.94 

TABLE  4b 


Amplification  Against  Wavelength  for  Mach  Number  0.5 


M  * 

oo 

0.5 

expCaCj-) 

X  =  2tt/oi 

1.0000 

6.54 

1.0164 

6.93 

1.0304 

7.34 

1.0482 

8.01 

1.0653 

8.94 

1.0779 

10.05 

1.0845 

11.00 

1.0891 

12.15 

1.0909 

13.01 

1.0916 

13.63 

1.0919* 

14.65 

1.0916 

15.41 

1.0911 

16.22 

1.0901 

17.11 

1.0889 

18.08 

1.0873 

19.13 

1.0855 

20.30 

1.0833 

21.59 

1.0809 

23.01 

1.0783 

24.60 

1.0753 

26.38 

1.0722 

28.39 
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TABLE  Ac 


Amplification  Against  Wavelength  for  Mach  Number  1.0 


H  = 

oo 

1.0  j 

exp(aCT) 

1  i 

X  =  271 /a  1 

I 

1.0000  j 

7.39  i 

1.0129 

7.81  j 

1.0267 

i 

8.36  j 

1.0361 

8.85  j 

1.0445 

9.38 

1.0517 

9.97 

1.0597 

10.84 

1.0657 

11.85 

1.0700 

13.03 

1.0723 

14.40 

1.0731* 

16.03 

1.0729 

17.96 

1.0724 

17.98  ! 

1.0715 

1 

19.12  1 

1.0702 

1  20.37  | 

1.0686 

21.77  , 

1.0667 

!  23.34  j 

1.0645 

|  25.12  j 

1.0620 

27.14 

1.0593 

j  29.46 

1.0562 

!  32.15 

1 

i 

1 

TABLE  Ad 


Amplification  Against  Wavelength  for  Mach  Number  2.0 


-0*1 


FIGURE  t.  PLOT  OF  v(y)  =  0-5 ( 1  ♦  TANH(y)). 


FIGURE  2.  THE  VARIATION  OF  THE  NEUTRAL  POINT  PLOTTED  AS  a  AGAINST  (1)  CR  AND  (2)  M 


=  AMPLIFICATION  FACTOR 


FIGURE  3.  AMPLIFICATION  AS  A  FUNCTION  OF  WAVELENGTH  FOR  MACH  NUMBERS 

0,0-5, 1-0  AND  2-0. 


PLOTTED  FOR  MACH  NUMBERS  0-0  >0,1*5  AND  2*0 


FIGURE  A 


IMAGINARY  PART  OF  THE  FUNCTION 


FIGURE  5 
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